Spontaneous and triggered vortices in polariton OPO superfluids 
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We study non-equilibrium polariton superfiuids in the optical parametric oscillator (OPO) regime 
using a two-component Gross-Pitaevskii equation with pumping and decay. We identify a regime 
above OPO threshold, where the system undergoes spontaneous symmetry breaking and is un- 
stable towards vortex formation without any driving rotation. Stable vortex solutions differ from 
metastable ones; the latter can persist in OPO superfluids but can only be triggered externally. 
Both spontaneous and triggered vortices are characterised by a generalised healing length, specified 
by the OPO parameters only. 
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Since the first observation of stimulated scattering [l| , 
resonantly driven polariton microcavities have been the 
subject of intensive research. Significant advances have 
taken place towards a new generation of low threshold 
lasers and ultrafast optical amplifiers and switches. How- 
ever, only very recently resonantly pumped polaritons 
have been shown to exhibit a new form of non-equilibrium 
superfluid behaviour In the OPO regime (Hj, po- 

laritons are continuously injected into the pump state 
and, above a pump strength threshold, undergo coher- 
ent stimulated scattering into the signal and idler states. 
Superfluidity has been tested through frictionless flow 
by triggering a traveling signal with an additional pulsed 
probe laser. Moreover, metastability of quantum vortices 
and persistence of currents have been proven by using a 
pulsed Laguerre-Gauss beam Vorticity has been ob- 
served to be transferred into to the OPO signal and to 
persist in absence of the driving rotating probe. 

The polaritonic system is intrinsically non-equilibrium: 
Continuous pumping is needed to balance the fast polari- 
ton decay, of the order of picoseconds, and maintain a 
steady state regime. In strong contrast with equilibrium 
superfluids, which ground state is flowless, pump and 
decay lead to currents that carry polaritons from gain 
to loss dominated regions. As a result of the interplay 
between these currents and a confining potential, polari- 
tons non-resonantly injected into a microcavity have been 
shown to become unstable to spontaneous formation of 
vortices d, Q and vortex lattices 0]. For resonant ex- 
citation, currents arise in the OPO regime due to the 
simultaneous presence of pump, signal and idler emit- 
ting at different momenta (see Fig. [T|). In this Letter we 
identify a regime of the OPO, where remarkably, even in 
the absence of disorder or trapping potentials, the system 
undergoes spontaneous breaking of the system symmetry 
and becomes unstable towards the formation of a quan- 
tized vortex state with charge m ~ ±1. We show that 
these spontaneous stable vortex solutions are robust to 



noise and to any other external perturbation, and thus 
should be experimentally observable. Spontaneous sta- 
ble vortices differ from metastable ones, which can only 
be injected externally into an otherwise stable symmetric 
state, and which persistence is due to the OPO superfiuid 
properties 01 ■ The metastable vortex is a possible but 
not unique stable configuration of the system. We find 
that shape and size of the metastable vortices are inde- 
pendent on the external probe. In addition, like in equi- 
librium superfluids, both stable and metastable vortices 
are characterised by a healing length which is determined 
by the parameters of the OPO system alone. Controlled 
creation of vortices in OPO has been recently achieved 
by a weak continuous probe Q . Metastable vortices have 
been also recently discussed for non-resonantly pumped 
polariton condensates in Ref. [l^l ■ 

Model We describe the OPO dynamics via Gross- 
Pitaevskii (GP) equations (Tlj for coupled cavity and cx- 
citon fields '0cx(r, i) with pumping and decay(?i = 1): 
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We neglect the exciton dispersion and assume a quadratic 



dispersion for photons, ujc 



O 2mc 



. The fields decay 



with rates kx,c and flft is the Rabi splitting. The cavity 
field is driven by a continuous wave pump. 



Fp(r,i) = J'/^,^^(r)e' 



;(kp-r— tiipt) 



(2) 



where J'fj,,(Tp is either a Gaussian or a top-hat spatial 
profile with strength fp and full width at half maximum 
(FWHM) dp. The exciton interaction strength gx can 
be set to one by rescaling both fields 4'x,c and pump 
strength Fp by \/^r/ {2gx)- For the simulations shown 
in this Letter mc = 2 x lO^^mo, the energy zero is fixed 



2 



C \ \ 

\ 

\ 
\ 

_.x.... : 








idler _ 








si 


LP 





-2-101234 



FIG. 1: (Color online) OPO spectrum for a top-hat pump of 
FWHM ap = TO^m and intensity /p = 1.24/^*''^ where /^""^ 
is the threshold strength for OPO. Polaritons injected reso- 
nantly at (fcp,0) and lo^ undergo coherent stimulated scatter- 
ing into the signal and idler states, which are blue-shifted with 
respect to the bare lower polariton (LP) dispersion (green 
dotted line) because of interactions. Cavity photon (C) and 
exciton (X) dispersions are plotted as gray dotted lines. 
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to ujx = i^c- (case of zero detuning), f2/^ ~ 4.4mcV and 
*^x.c ai'C fixed so that to give a polariton lifetime of 3ps. 

For homogeneous pumps, the conditions under which 
a stable OPO switches on can be found analyticall y by 
making use of the plane- wave approximation [l2l . Il3| . 
However, for a finite si ze p ump ([2|). one has to resort 
to a numerical analysis [ll] . Here, we numerically solve 
Eq. ([1]) on a 2D grid by using both a 5*^-order adaptive- 
step Runge-Kutta algorithm and the Crank-Nicholson 
method. Fixing the pump momentum (fcp, 0) close to the 
LP inflection point, we find the pump strength threshold 
/p*^ for which signal and idler states get exponentially 
populated. The OPO signal spatial profile |e"^cr,x 
can be obtained by either filtering in a cone around the 
signal momentum at a given time or by filtering the time 
resolved spectrum in an interval around the signal en- 
ergy (see Fig. [1]). The photon component of the filtered 
OPO signal profile and its supercurrcnt Vi^^ are 
shown in the first left panel of Fig. [2] We only select 
OPO solutions which reach a steady state. Note that 
the pump direction (fcp,0) leaves the symmetry y i— ^ —y 
intact. This symmetry, while allowing vortex-antivortex 
pairs, does not permit single vortices. 

One can perform a dynamical stability analysis of the 
OPO by adding small fluctuations to the steady state 
mean-field: The existence of modes with positive imag- 
inary part in the excitation spectrum indicates dynam- 
ical instability towards the growth of different modes. 
Equivalently, stability can be checked numerically by in- 
troducing a weak noise. In particular, we add white 
noise as a quick pulse to both modulus |-!/'x,c(k, t)| and 
phase (/)x,c(k, i) of excitonic and photonic wavefunctions 
in momentum space. The noise added to the phase has 
amplitude 27r, and for the modulus we specify the noise 



FIG. 2: (Color online) Filtered signal profile |7/)J(r, f)| with 
supercurrents V(j!>p (r, t) (left panels) and full momentum 
emission IV'cCi^a;, 0, f)| (right in arb. units) at three differ- 
ent times: t = ^ (first row), t — 432ps (second) and 1.2ns 
(third). The strength of the cTp = 35^m top-hat pump is 
fp — 1.12/p"^'. At t = a pulsed weak random noise of 
strength 0.01 (see text) is added to the OPO steady state 
(first row) and a,t t — 432ps a vortex, with m — —1, enters 
the signal and settles into a steady state. The noise strength 
only influences the time required for the vortex to appear. 



strength in units of the maximum of the pump intensity 
in momentum space. 

Stable vortex solutions Remarkably, we have singled 
out steady symmetric OPO states, as shown in the first 
row of Fig. [21 which are unstable towards the spontaneous 
formation of stable vortex solutions. Once the y i— > —y 
symmetry is broken by the noise pulse of any strength, we 
have observed a vortex with quantised charge to = ±1 
(to = =f1) entering and stabilising into the OPO sig- 
nal (idler). In the case of Fig. [2] and the right panel of 
Fig. El the noise strength is 0.01 and 432ps after the noise 
pulse, a vortex with to = —1 (m = -1-1) enters the signal 
(idler) and stabilises. Different noise strengths do not 
affect the final steady state, but only the transient time 
the system needs to accommodate the vortex and reach 
the steady configuration. Note that parametric scatter- 
ing constrains the phases of pump, signal and idler by 
20p = (ps ~\~ 4>i- Therefore a vortex in the signal implies 
an antivortex in the idler and viceversa. 

Further, we examine whether this vortex steady state is 
dynamically stable by applying an additional noise pulse. 
For weak noise, with a strength up to 0.1, the vortex is 
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FIG. 3: (Color online) Steady state filtered signal profile (dot- 
ted line) ipQ{x,0,t) for y ~ before the arrival of either a 
vortex probe Q with apt — l/^m (left panel, red dashed line) 
or a noise pulse of strength 0.01 (right panel, red dashed line) 
— same OPO conditions as Fig. [51 After the arrival of any 
perturbation breaking the y ^ —y symmetry, the same vor- 
tex with charge m = ±1 (solid shaded curve) stabilises into 
the signal. 

stable and can only drift around a little before settling 
again into the same state. For strong noise, with strength 
1 and above, the vortex gets washed away, but after a 
transient period, the very same state enters and stabilises 
again into the signal, with the possibility of flipping vor- 
ticity, as discussed later. Different noise strengths do not 
affect the final steady state, but only the transient time. 

Alternatively, one can break the y t-^ ~y symmetry by 
a weak pulsed vortex probe, and assess whether the stable 
steady state is in any way dependent on the external 
perturbation. As already shown in Refs. 0, [l3l, vortices 
with charge m = ±1 can be generated in the OPO signal 
and idler, by adding a Laguerre-Gauss pulsed probe: 

where the probe momentum kph and energy Wpf, are reso- 
nant with either the OPO signal or idler state. The phase 
ip winds from to 2tt around the vortex core Vpi,. Shortly 
after the arrival of the pulsed probe (o-j = Ips) there 
are two possible scenarios: either the vortex is imprinted 
into the signal and idler and drifts around, or no vortex 
gets transferred. However, the homogeneous OPO states 
which are unstable towards the spontaneous formation of 
stable vortices following a white noise pulse, exhibit the 
same instability following a vortex probe pulse (see left 
panel of Fig. [3]). The steady state vortex is independent 
on both the probe intensity fph and size apb, however the 
weaker the probe the longer the vortex takes to stabilise, 
between 30 and 400ps for our system parameters. The 
stable vortex following the Laguerre-Gauss probe is ex- 
actly the same as the one triggered by a weak white noise 
(see right panel of Fig. indicating that the probe acts 
only as a symmetry breaking perturbation. 
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FIG. 4: (Color online) Generation of an m = -fl metastable 
vortex solution into the OPO signal {fp — 1.24/p'' and 
(7p = 70/^m). First row: OPO filtered signal spatial profile 
together with currents (left) and interference fringes (right) 
are plotted at t = — 24ps before the arrival of either a strong 
enough vortex probe Q or a strong enough noise pulse. The 
metastable vortex lasts for as long as our simulation (last 
row, t = 1.8ns) and requires a threshold in the intensity of 
the perturbation breaking the y ^ —y symmetry. The vortex 
appears in the interference fringes as a fork-like dislocation. 



We can therefore infer that there are OPO conditions, 
as the one of Fig. [2l where the y i— >■ —y symmetric solu- 
tion is dynamically unstable, and any symmetry break- 
ing perturbation allows the signal and idler to relax into 
a stable steady state carrying a vortex with charge ±1. 
This suggests that such a vortex is the genuine unique 
OPO stable state, which however cannot be accessed 
without breaking the y i— > —y symmetry. Instability of 
the uniform state to spontaneous pattern (e.g., vortex) 
formation is a typical feature of systems driven away from 
equilibrium [l5|. Similarly we find conditions for which 
the uniform OPO solution is unstable to spontaneous for- 
mation of a quantised vortex. 

The system symmetry can also be broken by spatial 
disorder. Indeed, as for non-resonantly pumped polari- 
tons 0, we have found in our OPO simulations sponta- 
neous stable vortices forming in a disordered landscape. 

Metastability In addition to stable vortices, we have 
found OPO conditions supporting metastable vortex so- 
lutions. In this case, the symmetric OPO steady state is 
dynamically stable, but, because of its superfiuid proper- 
ties, can support persistent metastable currents injected 
externally. Metastable solutions can be equally induced 
by cither a vortex probe ([3]) or a noise pulse. However, 
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differently from the stable case, such a solutions require 
a threshold in the perturbation breaking the system sym- 
metry. One example is shown in Fig. 2] the steady state 
shown in the second row can be triggered by adding a 
white noise pulse to the symmetric OPO steady state 
(first row of Fig. 0]) only for noise strengths larger than 
0.1. For weaker noises, the OPO signal is slightly per- 
turbed and rapidly goes back to its steady state vortex- 
less configuration. Alternatively, a vortex probe ([3]) can 
also be applied. For example, for a probe vortex with 
Tph ~ (— 6,— 5)/xm, (Tpf, = 4.5/im, which is resonant with 
the OPO signal momentum and frequency, a probe in- 
tensity /pb > 0.45/p is required in order to generate the 
steady state vortex shown in the last row of the figure. 
Right panels of Fig. |4] show the interference fringes be- 
tween signal and pump state, obtained by considering the 
full emission in space |^/'c(r, t)\. A vortex corresponds to 
a fork-like dislocation in the interferences. 

The spatial position of both stable and metastable 
steady state vortices is close to the position where the 
OPO signal has the currents pointing inwards (see Figs. [2] 
and 2]) . We find stable vortices appearing quite close to 
the pump threshold. At higher powers, the OPO tends 
to switch off in the middle of the excitation spot and the 
system tends not to reach a steady state. Finally, we 
checked that m = ±1 (m = =f1) vortex solutions can 
appear only into the OPO signal (idler). A vortex probe 
pulse of any charge m injected resonantly to the pump 
momentum and energy gets immediately transferred to 
an m = ±1 (rn = =f1) vortex in the signal (idler), leav- 
ing the pump vortexless. The stability of signal m > 1 
vortices is outside the scope of this work [3l ■ 

Conservation of charge When generated by a noise 
pulse, stable and metastable vortices have equal proba- 
bility to be triggered with either charge ±1. Similarly, 
when vortices are triggered via a vortex probe, their vor- 
ticity can flip during the transient period. Often, the 
flipping is caused by the spontaneous appearance of two 
antivortices at the edge of the signal, one recombining 
with the triggered vortex. Topological charge inversion 
has been already predicted to occur in confined atomic 
Bose-Einstein condensates at intermediate interactions 
strengths [l7| : In the presence of an asymmetry, breaking 
rotation invariance, vorticity is not a conserved quantity. 

Healing length An approximate analytical expression 
for the vortex healing length can be derived for homoge- 
neous pumping, assuming that only signal and idler carry 
opposite angular momentum m, 

while the pump remains in a plane-wave state, ipP{r) = 
y^gjkp-r^ as supported by our numerical analysis. For 
pump powers close to OPO threshold, it can be shown 
that signal and idler steady state spatial profiles are 
locked together and satisfy the following complex GP 



equation 

2mc' \ dr'^ r dr r-^ J ' 

where jaj ~ ffx Keeping aside that a is a complex 

parameter, this equation describes a vortex profile [isj 
with healing length 

^ = (2mc<?x V^)''^' • (4) 

This expression is similar to the one of an equilibrium 
supcrfluid, with the condensate density replaced by the 
square root of the product of the signal and idler densi- 
ties. Further above threshold, the signal and idler profiles 
are no longer locked together, and they start to develop 
different radii. In both the simulations of Figs. [3] and SI 
we find 5 ~ 4^m, compatible with the estimate (|4]). Re- 
cently, the controlled creation of OPO vortices by a weak 
continuous probe Q has allowed to experimentally test 
the validity of Eq. (g]). 

To conclude, we have shown that close to the OPO 
threshold, the polariton supcrfluid can spontaneously 
break the y i— >■ —y system symmetry and, even in ab- 
sence of driving rotation, trapping, or disorder potential, 
can be unstable towards the creation of a quantised vor- 
tex state of charge ±1. While the OPO symmetric state 
is generally stable, metastable vortices can be injected by 
a strong enough external probe. The size of both types 
of vortices is given by a healing length, which, close to 
threshold, is analogous to the one of equilibrium super- 
fluids. Long lived metastable vortices have been exper- 
imentally realised in Ref. Q, providing an evidence for 
existence of persistent currents in this system. Since the 
parameters of our simulations are close to those of cur- 
rent semiconductor microcavities, the existence of stable 
spontaneous vortices in OPO should also be within exper- 
imental reach. Moreover, these stable vortex states are 
very robust to noise. Thus, apart from being yet another 
example of exotic behaviour in polariton superfluids, they 
carry potential for applications in, e.g., quantum infor- 
mation. 
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